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Abstract 

We study conductance spectroscopy of a two-dimensional junction between a normal metal and 
a strongly- correlated superconductor in an applied magnetic field in the Pauli limit. Depending on 
the field strength the superconductor is either in the Bardeen-Cooper-Schrieffer (BCS), or in the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state of the Fulde-Ferrell (FF) type. The strong correla- 
tions are accounted for by means of the Gutzwiller method what leads naturally to the emergence 
of the spin-dependent masses (SDM) of quasiparticles when the system is spin-polarized. The case 
without strong correlations (with the spin-independent masses, SIM) is analyzed for comparison. 
We consider both the s-wave and the d-wave symmetries of the superconducting gap and concen- 
trate on the parallel orientation of the Cooper pair momentum Q with respect to the junction 
interface. The junction conductance is presented for selected barrier strengths (i.e., in the contact, 
intermediate, and tunneling limits). The conductance spectra in the cases with and without strong 
correlations differ essentially. Our analysis provides thus an experimentally accessible test for the 
presence of strong- correlations in the superconducting state. Namely, correlations alter the distance 
between the conductance peaks (or related conductance features) for carriers with spin-up and spin- 
down. In the uncorrelated case, this distance is twice the Zeeman energy. In the correlated case, 
the corresponding distance is about 30-50% smaller, but other models may provide even stronger 
difference, depending on details of the system electronic structure. It turns out that the strong 
correlations manifest themselves most clearly in the case of the junction with the BCS, rather than 
the FFLO superconductor, what should make the experimental verification of the present results 
simpler. 

PACS numbers: 74.45. +c, 71.27.+a, 71.10.Ca, 74.50. +r 
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I. INTRODUCTION 



The search for evidence of strong electron correlations in condensed matter has con- 
centrated in recent years on superconducting state in unconventional materials and its co- 
existence with magnetism. One of such examples is the search for experimental evidence 
for the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state. The FFLO state 
was proposed theoretically in the IDGOs.-"^ In this unconventional superconducting state the 
Fermi wave vector difference for the electrons with spin-up and -down due to the presence 
of Zeeman term makes it favorable for the Cooper pair to acquire a nonzero total momen- 
tum Q = 2q. Consequently, the phase of the superconducting gap parameter oscillates 
spatially with the wave vector Q. By forming such a condensate of moving Cooper pairs, 
the superconducting state persists to magnetic fields remarkably higher than the Pauli Hc2 
limit. The FFLO state has suddenly gained renewed interest recently (for a review see 
Ref. [sl) because of its possible detection in the heavy fermion superconductor CeCoIns,-^- 
although the nature of the high-field low-temperature phase observed in this system is still 
under an intensive debate after antiferromagnetism has been observed in the vicinity of 
this phase.-^^-i^ The FFLO state has also been proposed for k-(BEDT-TTF)2Cu(NCS)2,— 



Ref. 



T)2SF5CH2CF2S03,— and other layered organic superconductors (see References in 



3). 



Also, existence of the FFLO state has been indicated in other heavy-fermion sys- 



tems: PuRhGas,— Ce2PdIng^ (see Ref. Il8|, Sec. V.B.I for a more detailed account), as well 
as in the pnictide superconductor LiFeAs.— The FFLO state has also been investigated in 
high density quark and nuclear matter,— as well as in optical lattices.—"— 

All the systems considered so far to be a host to the FFLO phase have a reduced di- 
mensionality, what is crucial for the FFLO phase stability, as then the orbital effects are 
suppressed and the Pauli effect (Zeeman splitting) may become the dominant factor. An- 
other obvious feature, which suppresses the orbital effects, is the heavy quasiparticle mass. 
These characteristics of possible FFLO hosts indicate that these systems are likely to have 
strong electron (fermion) correlations and thus also possess specific features resulting from 
them. 

The role of strong correlations in the most likely candidate for the FFLO state, CeCoIns 
is essential not only because this system is a heavy fermion superconductor, with very 
narrow bands originating from 4/ electrons hybridized with bd — 6s states. What is equally 



important, the spin-dependent effective masses (SDM) of quasiparticles have been directly 
observed in this systent^ by means of the de Haas- van Alphen oscillations in a strong applied 
magnetic field. SDM are one of the hallmarks of strong correlations, as they appear naturally 
in theories incorporating correlations (Gutzwiller,— slave-bosons,—!^ dynamical mean field 
theory,— fluctuation-exchange approximation^^) , when the system is spin-polarized.— 

Because of the above reasons, it is important to study the effect of correlations on the 
FFLO phase. Such analysis has already been performed in a few cases,—""— and it indicates, 
among others, that the interelectronic correlations play an important role in forming and 
stabilizing the FFLO phase. 

In the present paper we concentrate on providing an experimentally-accessible concrete 
characteristics of a superconducting state with strong correlations. Namely, we study con- 
ductance of a normal metal - superconductor junction (NSJ) with the strongly-correlated 
superconductor in either the Fulde-Ferrell (FF) type of the FFLO state, or the Bardeen- 
Cooper-Schrieffer (BCS) state (the latter is stable in lower fields). Conductance spectroscopy 
of such junction is an experiment sensitive to both the phase and amplitude modulation of 
the superconducting order parameter, and therefore it is a candidate technique for providing 
a direct evidence for the presence of the FFLO phase. In that situation, a crucial role is 
played by the Andreev reflection (AR) processes.— In the simplest view of the Andreev re- 
flection, an incident electron entering from the normal metal into the superconductor (SC) is 
converted at the NSJ interface into a hole moving in the opposite direction (to the incident 
particle) and Cooper pair inside SC. Such processes increase conductance of the junction (in 
an ideal case by a factor of two), which is analyzed in the framework provided by Blonder, 
Tinkham, and Klapwijk.— 

The conductance characteristics for a NSJ with superconductor in the FFLO state has 
already been investigated for both the cases of the FF (with A(r) = Aqc'^"")^""— and the 
Larkin-Ovchinnikov (A(r) = Aq cos(Qr))^ types of FFLO states, as well as for the case of 
superconductor with a supercurrent^'^ (i.e. the situation similar to ours from formal point 



of view). See also Refs. 



44 



- 471 for the case of NSJ with BCS state of the d-wave symmetry. 
None of the above papers have taken into account strong electron correlations. 

Here we consider both the cases of s-wave and rf-wave strongly-correlated superconduc- 
tor in magnetic field and in the Pauli limiting situation (i.e., we neglect orbital effects, as 
the Maki parameter— in the systems of interest is quite high^) . The strong correlations are 



taken into account by assuming dispersion relations with SDM of quasiparticles and with the 
correlation field, as given e.g. by the Gutzwiller approximation,— or slave-boson theory.— 
The case without strong correlations (with spin-independent masses, SIM) is analyzed for 
comparison. In low magnetic fields the superconductor is in the BCS state, and in higher 
magnetic fields a transition to the FFLO state takes place. We consider only the simpler 
FF type of FFLO state as we intend to single-out novel features of the situation with strong 
correlations in the simplest case (the analysis of LO state is much more complex—). Our 
study already leads to interesting, novel results. We set the direction of the Cooper pair 
momentum Q as either perpendicular, or parallel to the junction interface, with more atten- 
tion paid to the latter situation. The analysis is performed in a fully self-consistent manner. 
Namely, we select Cooper pair momentum Q minimizing the free energy of the system and 
we determine the chemical potential /i in each phase separately so that the particle num- 
ber n is kept constant. Such an adjustment of is required even for the BCS state for 
the narrow-band case. Also, such a careful examination of the superconductor properties 
is important, and non-self-consistent calculations may lead to important alterations of the 
conductance spectrum.— 

As we deal with heavy quasiparticles on the superconducting side of NSJ, we should in 
principle take into account the Fermi-velocity-mismatch effects. Under those circumstances, 
the AR processes would be severely limited by a high effective barrier strength Z. On the 
other hand, AR is clearly observed in junctions with heavy-fermion superconductors^"^ and 
theoretical efforts have been made to understand why this is the case.—"— Based on these 
studies, we disregard the Fermi-velocity mismatch by assuming equal chemical potentials 
and equal average masses of quasiparticles on both sides of the junction. Namely, we choose 
masses on the normal side as mav, and on the superconductor side we have that (m-|- -|- 
m^)/2 = mavi with = 100 mo (where mo is the electron mass in vacuum), which 
roughly corresponds to the heaviest band of CeCoIns.— This assumption is, in our view, a 
justifiable simplification, as we would like to single out the novel features in their clearest 
form. Note also, that we consider a model situation with its parameters taken from the 
experiment for CeCoIns. 

In brief, we study conductance of NSJ with superconductor exhibiting strong electron 
correlations (SDM case). To single out novel features of such situation, we also study the 
uncorrelated case (SIM) and compare those results. 
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The paper is organized as follows. In Sec. [TTl we discuss the superconducting state of 
quasiparticles with SDM and SIM for a two-dimensional electron gas. In Sec. Illll we present 
the theory concerning conductance of a normal metal - strongly-correlated superconductor 
junction. In Sec. II VI we show conductance spectra for the cases with SDM and SIM. In Sec. 
|V]we discuss relation of our results to experiments and suggest their possible experimental 
verification. Finally, in Sec. |VT]we provide a brief summary. 



II. FULDE-FERRELL SUPERCONDUCTING STATE BASIC CHARACTERIS- 
TICS: MODEL AND METHOD 



As said above, here we consider a two-dimensional system of paired quasiparticles in 
the situations with SDM and SIM. The system of self-consistent equations describing such 
superconducting state has already been presented in detail in Refs. |3l| and |32l For the sake 
of completeness, we provide here a brief summary of our procedure. We start with the 
Hamiltonian 



"H — ^ ^kCToL«k(T + ^ ^k,k'Ok+q|Ot-k+q4'^-k'+q4a-k'+qt + —mKor, (l) 
kcr kk'q 

where Q = 2q is the wave vector of the Cooper pair center of mass, n = + ni is the band 
filling, m = n-|- — is the spin-polarization of the system, and is the total number of 
particles. The dispersion relation for the cases with SDM and SIM is chosen, respectively, 

as 



ko- 



JSIM) 

Skcr 



2mcr 



a{h + he 
- ah — fi, 



(2) 

(3) 



where h = gfiBH/2 with H being the applied magnetic field. The quantity hcor is the 
correlation field which appears naturally in both the slave-boson theory (it is equivalent 



to -(3 of Ref. 



261 1 and Gutzwiller approximation if this approximation is performed with 



care.— — Justification of a Hamiltonian with both the pairing part and SDM can be found 
(Appendix A) and in Ref. |57|. The spin- dependent quasiparticle mass is equal 



33 



in Ref. 

to rricr = niB/qcrin^rri), where is the bare band mass and q^{n,fn) is the band-narrowing 



factor. Explicitly (in the Hubbard f/ — )■ oo limit) the quasiparticle masses are given 
m„ 1 — rifj 1 — n/2 fn _ 1 

triB 1 — n 1 — n 2(1 — n) ms 



[rrir, 



cTAm/2), 



(4) 
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with Am = mi—m^. Next, as in the BCS theory, we take the pairing potential in a separable 
form and assume it is nonzero in a small region around the Fermi surface (for details see 



Refs. 



31 



32 



and 



581) 



Vk,k' = -K)?7k'7k', (5) 

where r/k = cos (akx) — cos (aky) for the d-wave case (with a = 4.62 A being the lattice 
constant for CeCoIn^^) and r/k = 1 for the s-wave case. Under such assumptions, the 
superconducting gap can be factorized as 

Ak,Q = AQ?7k. (6) 

Following the standard mean- field approach to Hamiltonian ([T]), we obtain the generalized 
free-energy functional J-" and the system of self-consistent equations as follows^i^^ 

T = -ksT V ln(l + e-^^-) + V(e[') - ^k) + N-^ + fiN + -m/i,,,,, (7) 
t k Vo n 



^ dfn N ^ dm V E] 

kcr k 



n 
N 



m = ^^a/(Ek.), (9) 



kcr 



- ^ Aq, (10) 



n 
N 



n = n^ + H = l Y.^ulf{E^^) + vl[l - /(i^k,-.)]}, (11) 

kcr 



where J^{T, H, /i; m, h^or, Aq, n) is the system free-energy functional for the case of a fixed 
number of particles^S (we fix the band filling at the value n = 0.97), Vq is the interaction 
potential, Uk, fk are the Bogolyubov coherence coefficients, f{E]icr) is the Fermi distribution, 
and Ucr is the spin-subband filling. The physical solution is that with a particular Q which 
minimizes the free energy F, which in turn is obtained from J-" by evaluating the latter at 
the values of parameters being solution to Eqs. (l8l)- (fTT!) . The state with Q = is called the 
BCS state, and that with Q 7^ - the FF state. 

The quasiparticle spectrum in the paired state is characterized by the energies (cf. also 
Ref.leO) 



E^. = E^ + ai\^\ Ek = V^t'^' + I Ak,Q|2, (12) 

^k ^ = 2 (^k+qt + ^-k+q;), ^k"^ = 2(^k+qt " ^-k+q^)- (13) 



Eqs. (l8i)- (fTT|) are solved by numerical integration over the reciprocal space. We use proce- 
dures from GNU Scientific Library^ as solvers. For the SIM case hcor = and we solve only 
Eqs. (pi)- fllll) . The numerical procedure has been elaborated in detail elsewhere.— Here, for 
completeness, we also provide in Tables 1 and 11 the numerical values of selected parameters 
for the situations with the s-wave and the d-wave symmetries of the superconducting gap, 
respectively. The quantity F^s is the free energy of the normal state, and therefore AF 
is the condensation energy. Also, Am. = m2 — mi is the mass difference and hcorFS is the 
correlation field value in the normal state. The free energies are calculated per elementary 
cell. The numerical accuracy is not smaller than on the level of the last digit specified. 



Table I. Equilibrium values of mean- field variables and related quantities 
for the s-wave solution with H = 10.01 T and T = 0.02 K. 


Variable 


Value 


Variable 


Value 


m 
hcor (K) 

Aq (k) 

/i(K) 
F(K) 


0.0129431 
-3.08230 
1.38922 
126.287 
61.18200288 


Am (mo) 

hcor FS (K) 

IQI(A"') 

\Q\/AkF 

AF {K)=Fns-F 


2.51322 
-3.26546 
0.00947 
1.08 
-0.00111351 



Table 11. Equilibrium values of mean-field variables and related quantities 


for the d-wave solution with H = 20.01 T and T = 0.1 K. 


Variable 


Value 


Variable 


Value 


m 


0.0268690 


Am (mo) 


5.21729 


hcor (K) 


-6.40870 


hcorFS (K) 


-6.53133 


Aq (K) 


1.27455 


IQI(A~^) 


0.0183 


/i(K) 


126.416 


\Q\/AkF 


1.15 


F(K) 


61.04342125 


AF (K) = Fns-F 


-0.00142436 



The input parameters in our method have the following values: the band filling n = 0.97, 
the lattice constant a = 4.62 A, the interaction potential strength Vo/n = 90 K (rf-wave) 
and Vo/n = 110 K (s-wave), the interaction potential width (cutoff) huc = 17 K, the 
quasiparticle average mass m^^ = 100 mo. The other parameters (in particular: m, hcor, 
Aq, fj,, Q, and 6'q) are determined from the solution procedure. 
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Phase Diagram: SDM, s-wave 
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Phase Diagram: SIM, s-wave 
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FIG. 1. (Color online) Phase diagram for the s-wave gap symmetry in the SDM (a) and SIM (b) 
cases. Light (yellow) region corresponds to Q = (BCS phase), the darker to the state with Q 7^ 
(FF phase) and the white to normal state (NS). Note the greater difference between SDM and SIM 
cases than for d-wave gap symmetry (see Fig. [2]) . 



Exemplary phase diagrams obtained on the applied field H and temperature T plane are 
exhibited in Figs. [Hand [2] for the s-wave and the d-wave cases, respectively. The angle is 
the angle between the maximum-gap (antinodal) direction and the Cooper pair momentum 
Q. Note that in both situations the FF state is more robust (i.e. the FF state fills a wider 
field-temperature range on the phase diagram) in the SDM case than in the SIM case. The 



mechanism of 
in Refs. 



31 



33 



;he FF state stabilization by strong correlations has been analyzed in detail 
For the sake of completeness, let us mention that this mechanism is based 
on a smaller Fermi wave vector splitting [Akp = kp'i^ — kp^) in the SDM situation. In such 
case, the system can resist more efficiently the destabilizing influence of the applied magnetic 
field (hence higher critical fields in the SDM case). Also, it turns out that the FF state can 
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a) 



Phase Diagram: SDM, d-wave 
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FIG. 2. (Color online) Phase diagram for the d-wave gap symmetry for the cases with SDM, (a) 
and SIM, (b). Light (yellow) region corresponds to Q = (BCS phase), the darker (blue, red) to 
the states with Q 7^ (FF phases) and the white to normal state (NS). The red regions correspond 
to the Cooper-pair momentum Q in the maximum-gap (antinodal) direction {0q = 0), whereas 
the blue ones to the momentum along the nodal direction (^q = vr/4). Note that this anisotropy 
results solely from the d-wave gap symmetry, as the unpaired gas is isotropic. The dashed line 
marks the BCS critical field Hc2 in the Pauli limit, and the dot-dashed line marks Hc2 for the 
solution with = 0. 



benefit to a greater extent than the BCS state from the smaller Akp, as the FF state has 
higher spin-polarization, which is necessary for the appearance of SDM (for details see Refs. 



31 



32, and 



58 



For further analysis of the Andreev reflection, we take the parameters obtained along 
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the T = 0.02K ^ line in Figs. [T] and |2l Therefore, the results will have, strictly 
speaking, practical relevance for T ^ Tgc, with the superconducting transition tempera- 
ture ~ 2 — 3 K, as can be seen from Figs. [T]and[2J 



III. JUNCTION CONDUCTANCE: THEORETICAL ANALYSIS 

For the analysis of the NSJ conductance we take the superconducting state parameters 
obtained self- consistently (from the procedure presented above). We consider here only two- 
dimensional NSJ for simplicity. Kinematics of the reflection may be analyzed by means of 
the Bogolyubov-de Gennes (BdG) equations^ 

Eu^{x) = -HouA^) + j rfx'A(s,r)tv(x'), (14) 
Etv(x) = -•Hot^a(x) + J rfx'A*(s,r)iz,(x'), (15) 

where s = x — x', r = (x + x')/2, and cr = ±1 is the spin quantum number of the incoming 
quasiparticle and Uo-(x) and fo-(x) are the particle and hole wave- function components. The 
one-particle Hamiltonian is given by 

^o(r) = -V-^V -ah- ah,or{v) - fi + \/(r), (16) 

where we have used the effective mass approximation^'^ to express the kinetic part as 



2m{x) 



V with m(r) = m{x) = mavQ{—x) + m„Q{x), similarly as in Refs. 



64 



67. The 



correlation field is nonzero only on the superconducting side of the junction (/icor(r) = 
hcor^ix)). Also, r = (x, y) and the interface scattering potential is chosen as a delta function 
of strength if, i.e. V{r) = H 6{x). The gap function can be Fourier transformed as following 

A(s,r) = j dke'^^A{k,r) = J dke^^^Ai^^Qe'^'' e{x), (17) 

with Ak,Q as in Eq. (|6]) but with the original set of coordinates rotated by a (cf. Fig. |3]). 
Explicitly, the superconducting gap we use from now on has the form (in the new coordinates) 

Ak,Q = Aq ^ cos {akx cos a — aky sin a) — cos {aky cos a + akx sin a)^ . (18) 

We neglect the proximity effects by assuming a step-like gap function. To solve the BdG 
equations we make the plane-wave ansatz. Namely, we assume that the two-component pair 
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wave function has the form 



^(r,a)= =6-^ ' , (19) 

v^{r)\w) J y?}e-"i'"|a) 

with u and v as constants and with a = —a (we have also dropped the a indices of u and 

v). We also remind that q = Q/2. By substituting (fT7|) and ( IT9|) into BdG equations ([1] 

( ITSl) and after some algebra we obtain the following matrix equation 

-E + ^k+q,a Ak,Q \ I u\a) 
Alk,Q -E - a-q,^ / \ V \a) 
where unpaired quasiparticle energies ^ko- are given by Eq. ([2]) or ([3]). Eq. fl20l) gives the 
dispersion relations for quasiparticles and quasiholes in the superconductor 



0, (20) 



where ^k'^''^'' have been defined in Eq. fll3p . One may check that the above equation is in 
accordance with Eq. (fT2l) . as -Ek+ = -E'kt (quasiparticle) and E^^ = — -Eki (quasihole) for 
incoming particle with spin a =t, as well as -E'k+ = -E-kj, (quasiparticle) and E^^ = — -E-kf 
(quasihole) for incoming particle with spin a =|. This holds as long as A^j^ q = A^_q, 
which is true for any real k. 

As already mentioned, we study the FF type of the FFLO superconducting state, in 
which A(r) = Aqc*^*''" and set the direction of the Cooper pair momentum Q = 2q 
as either perpendicular (Q = {Q,0)), or parallel (Q = (0,(5)) to the junction interface. 
The perpendicular configuration (Q = (Q,0)) may lead to accumulating of charge at the 
NSJ interface due to normal and/or supercurrent present in the FF state. Therefore we pay 
principal attention to the parallel configuration. Parenthetically, the accumulation processes 
are very slow for the case of heavy quasiparticles. 

As we consider electron injected from the conductor side of the junction (junction geom- 
etry is presented in Fig. |3]), the corresponding wave functions can be expressed as (we have 
omitted the spin part for clarity) 

^<(r) = I ] j e^*^"- + ° j '''''' + ^ ' '''''''' ^^^^ 
i'-M-M "'^'"^V-^ + cf |e-f, (23) 

12 




a) 



b) 



Normal metal, m 



av 



Superconductor, m^ 




FIG. 3. (Color online) a) Junction geometry for incoming particle of spin a =\. Normal-state 
and superconductor regions are marked. Interface lies at the x = line. The superconducting 
gap is also presented: a is the angle between the axis and maximum-gap direction. Full circles 
mark quasiparticles and empty ones mark quasiholes. Momentum of each of them is marked with a 
boldface letter, and amplitude with an italic letter. Namely, incoming particle has the momentum k, 
and amplitude 1, reflected hole has p and a, reflected quasiparticle: k', 6, transmitted quasiparticle: 
k2 , c, and transmitted quasihole: k^, d. The angle of incidence is equal to 6 and to the angle 
of reflection but other angles (of reflection of quasihole and those of transmissions) may differ (cf. 
also Fig. HI). In (b) and (c) we show explicitly the two d-wave configurations of the superconducting 
gap for the FF phase studied in the following: (b) corresponds to Fig. [7] and (c) to Fig. [HI 

where ip<{'r) and describe wave function on the normal-metal and superconductor 

sides, respectively. The quasimomenta (for quasihole) and kj' (for quasiparticle) are 
solutions of Eq. (12T|) for a given incident energy E propagating in the positive x direction. 
From the translational symmetry of the junction along the y direction comes conservation 
of the y momentum component. Namely, ky = k'y = Py = k^y = /c^. All the wave vectors 
are presented in Fig. |H 

We use boundary conditions with the appropriate masses^ and the interface potential 
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a) 



Normal state 



BCS Superconductor 



b) 




Normal state 



FF Superconductor, SIM 




FF Superconductor, SDM, h ~ h„, 
E. =E 



FIG. 4. (Colour online). The junction geometry in reciprocal space. All relevant vectors are 
marked. It can be seen that only the incident and reflection angles are equal to 9. It can be 
anticipated at this point that changing 9 for BCS state does not lead to drastic changes in the 
transmission/reflection probabilities, whereas for the FF state the situation is quite different since 
Q 7^ induces anisotropy in the reciprocal space. The energy E value has been chosen as 10 K 
for all graphs except (b) "FF Superconductor" for which E = 0.01 K ~ (for E > 0.5 K there 
would be no = £^k- regions in this case). The dashed lines are guide to eye and illustrate the 
conservation of momentum y-component. 



jump H] they cirG clS follows 



^/'<(r)|^=o = ^/'>(r)U=o, 



1 



1 di)^^ 



dx 



\x=0 



dx 



\x=0 



2H 



V'<(r)U=o- 



(24) 
(25) 



Those conditions lead to the following set of 4 equations^ for the amplitudes (a, 6, c, d) 
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l + b-cu2-dui= 0, (26) 

a — CV2 — dvi = 0, (27) 

iK{l-b) cu2i{q^ + kt^) duii{q^ + kt^) ^ 2H \ b) - (28) 
TUav m„ m„ It? 



which are similar to those in e.g. Ref. 



64L except in our case vectors are replaced by their 
X- components: e.g. k ^ k^, p ^ Px, and also SDM are properly accounted for (obviously 
in the SIM case we have that = = m^^). From the solution of Eqs. fl26l) - fl29|) one 
can obtain probabilities of the hole reflection p'^./^ = |ap^^^, particle reflection = 
quasiparticle transmission 



■|«2P \v2f 



pI = \c\^ma,^ - - (30) 

kx 

and quasihole transmission 

PtH = \d\'ma.-^ '"-^ ^ , ^""^ (31) 



where the a superscript indicates the spin of the incoming electron. In the following we use 
the dimensionless barrier strength Z = 2mavH /{kpfi^), where we define Fermi wave vector 
kp using the zero-field value kp = j^y/2mavfi- Note also that we do not use the assumption 
k = k' = p = kf = k2 ^ kp utilized at this point in majority of the papers on Andreev 
reflection spectroscopy, because we deal with heavy quasiparticles for which fi is of the order 
of 100 K. Therefore the usual assumption p ^ E is not, strictly speaking, applicable in the 
present situation. 

IV. RESULTS AND PHYSICAL DISCUSSION 

Differential conductance (G = dl/dV) can be obtained from the reflection and transmis- 
sion probabilities^!^ in a straightforward manner 

Gns = I rf^cos^il -pUE,e)+p-,,{E,e)]. (32) 

The final result of our calculation is the total conductance G averaged over spin and nor- 
malized with respect to the conductance G^^ of the junction with A = but still with the 
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same other parameters {m„, /x, hcor), as the superconducting state. Namely, 



+G^ 

G = ."^ . (33) 



This quantity is exhibited in the following figures, sometimes with the spin-resolved con- 
ductance = G'^g/G'^^. We assume the barrier strength equal to Z = (contact limit), 
Z = 0.5 (intermediate limit), and Z = 5 (tunneling limit). The case of Z = 5 refiects not 
only the situation for planar NSJ with a thick insulating layer, but also that encountered in 
Scanning Tunneling Spectroscopy (STS) experiments.— 

Our goal in the following is to identify novel, model-independent features of the strongly- 
correlated situation (i.e., with SDM). Namely, those features should not depend on the 
assumed dispersion relation or the pairing-potential strength. 



A. s-wave pairing symmetry 

In Fig. |5] the conductance for the s-wave gap symmetry and Q vector oriented perpen- 
dicular to the junction, is presented. It can be seen that there are peaks in the conductance 
originating from AR processes of quasiparticles having different spins, that take place when 
the e nerg y E of the incoming electron fits into the so-called Andreev Window (AW), see 



Refs. 



39 



(Fig. 3), 



40 



and Ref. 



58 



(Chapter 5, Figs. 5. Id, 5.4b) for more details. These 
peaks are separated by a distance equal to twice the Zeeman energy {2h = q^bH) only in 
the case without strong correlations (SIM). For the SDM case the correlations compensate 



the Zeeman splitting (by means of hcor and mo-, cf. Refs. |3ll and 1321 ) and as result the 
conductance peaks are closer than twice the Zeeman energy. We identify this feature as 
a hallmark of strong correlations in the superconducting state. Another interesting feature 
differentiating the SIM and SDM cases is absence of the a =t peak for SDM when mag- 
netic field if > 12 T. For such fields the junction is transparent to incoming particles with 
a =t because the Andreev window— falls below E = {). In other words, the quasiparticle 
energy £'k+ within FF superconductor is below zero around the whole Fermi surface. This 
leads to breaking of Cooper pairs and produces normal state region filling whole angular 
space around the Fermi surface (see Fig. Hb, SDM case). Since there are normal particles 
with a =t within the FF superconductor, the incoming a =t quasiparticle does not feel the 
superconducting gap presence, and the junction is transparent, what yields G-t- ~ 1. 
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b) FF, SDM, s-wave 




2 4 6 8 10 

Energy (K) 



FIG. 5. (Colour online). Conductance spectra for the case of s-wave FF state for SIM (a) and 
SDM (b) cases. The Q vector is oriented perpendicular to the junction and the intermediate 
barrier strength Z = 0.5 is taken. The value of the gap is (a) Aq = 0.11 K; (b) Aq = 0.36, 0.70, 
1.09, and 1.53 K for the decreasing magnetic field. The distance between the peaks is twice the 
Zeeman energy 2h = g^sH only for the SIM case. In the SDM case the correlations compensate 
the Zeeman splitting (by means of hcor and mo-), and the peaks are closer than gusH. 

In all the following figures the parallel orientation of the Q vector with respect to the NS J 
interface has been assumed. In Fig. [6] the NS J conductance for the s-wave gap symmetry 
has been presented. Again, at high magnetic fields if > 12 T the junction is transparent to 
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FIG. 6. (Colour online). Conductance spectra for the case of s-wave FF state for SIM (a-c) and 
SDM (d-f), the Q vector oriented parallel to the junction, and selected Z values. In (b) and (d) 
also the spin-resolved signals Ga are presented. The distance between the characteristic features 
is shown in (a) and (d). In (a) and (e) we provide the values of the gap Aq (they are identical 
in (a)-(c) and (d)-(f)). In the SDM case for H > 12 T there are no features of the spin- up signal 
because the junction is transparent for incoming electrons with spin a =t (for explanation see main 
text). 
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incoming quasiparticles with a =\. In the present case it is difficult to discern characteristic 
features of the conductance from the spin-up and spin-down channels in such a manner, that 
the splitting of the peaks could be measured. For this purpose, the spin-resolved signals 
would have to be singled out, as shown in Fig. Ebd, because the spin-specific features of 
the total conductance are subtle and could be smeared out at finite temperature or due to 
other effects (e.g. inelastic scattering). Again, the characteristic features of spin-up and 
spin-down signals are separated by a distance equal to twice the Zeeman energy for SIM 
(Fig. [6^) and are closer for SDM (Fig. |6li). 

B. d-wave pairing symmetry 

In Fig. [7] the conductance in the case of FF state with = is presented. Such phase is 
stable in the high- field regime (see Figs. [1] and |2]). Note that by fixing the direction of Q with 
respect to the NSJ interface we fix also the angle a (see Fig. [3]), as 6'q is determined from the 
results presented in Sec. [TTl Namely, the parallel vector Q orientation with respect to the 
junction interface implies that a = for 0q = (cf. Fig. (Hb) and a = 7r/4 for = 7r/4 (cf. 
Fig. Eb). In the case with = no remarkable, model-independent differences between 
the SDM and the SIM cases appear, as all peaks present in Fig. [7| come from a =J, electrons 
(see Fig. [7^, where the a =t signal has been plotted). 

The conductance spectra for the d-wave FF phase with 9q = n/A (with a = vr/4) have 
been presented in Fig. [HI As in the s-wave case, and for the same reasons, at high mag- 
netic fields the junction is transparent to spin- up quasiparticles in the SDM case. Only at 
H < 14.4 T we were able to discern characteristic, spin-specific features of the spectra (see 
Fig. [H^d for the spin- resolved spectra). These features are again split by twice the Zeeman 
energy for SIM, and are closer for SDM. To identify the spin-specific features, spin-resolved 
spectra have to be analyzed, similarly as in the s-wave case. 

Finally, in Fig. IHlwe show the conductance spectra for the d-wave BCS state with (100) 
contact. In this case, in the tunneling limit [Z = 5) the peaks originating from AR of 
quasiparticles with different spins, are most clearly visible. As previously, these peaks are 
split by twice the Zeeman energy for SIM, and are closer for SDM. We identify this case as 
the most promising for experimental verification, as discussed in the following. 
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FIG. 7. (Colour online). Conductance spectra for the case of d-wave FF state for selected Z 
values, for the SIM and SDM cases. The Cooper pair momentum is oriented along the maximum- 
gap (antinodal) direction (i.e. = and a = 0; cf. Figs. [2^ and [8)3). In (b) we provide the values 
of the gap Aq, and in (a) we plot also the spin-up conductance. The magnetic field is close to lic2- 
There is no clear distinct feature, which differentiates between the SIM and the SDM situations 
for this configuration. 

V. RELATION TO EXPERIMENT 

Our results imply that the splitting between the spin-up and the spin-down features of 
the conductance spectra is equal to twice the Zeeman energy only in the non-correlated case 
(SIM). In the strongly-correlated case, due to the presence of spin-dependent masses (SDM) 
vrifj and correlation field /icor, the separation of the spin- up and the spin-down features differs 
essentially. In the present case of a two-dimensional, correlated electron gas, this separation 
is smaller (because mo- and h^oT compensate the Zeeman term; typically h^or ~ 0.5 x (— /i). 



cf. Refs. 



31 



and 



321 ) . but in general it may be larger. For example, in the two-dimensional 
Hubbard model, our recent calculations^ yield typically h^or ~ 5 x /i, and therefore in that 
model correlations enhance splitting of the conductance peaks. 

It should be in principle possible to measure the conductance-peaks splitting experimen- 
tally. Especially, the BCS case with (100) contact and high barrier strength Z (Fig. [9bf) 
looks promising, as the peaks are clearly visible, and the BCS state exists in lower magnetic 
fields than FFLO, what should make the whole analysis simpler (the orbital effects,— which 
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FIG. 8. (Colour online). Conductance spectra for the d-wave FF state with 0q = 7r/4 (Q along 
the nodal direction, a = 7r/4; cf. Figs. [2^ and[3l:) for selected barrier strengths for the SIM (a-c) 
and the SDM (d-f) cases. In (b) and (e) we provide the values of the gap Aq. In (a) and (d) 
also the spin-resolved conductance has been presented to identify spectra features for both spin 
channels. These features are separated by twice the Zeeman energy for the SIM, and are closer 
again for the SDM case. 
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FIG. 9. (Colour online). Conductance spectra for the d-wave BCS state with (100) contact (a = 0) 
for selected barrier strengths for the SIM (a-c) and the SDM (d-f) cases. In (a) and (d) we provide 
the values of the gap Aq. In (c) and (f) in the tunneling regime {Z = 5) conductance peaks from 
spin- up and spin-down channels are clearly visible already in the total conductance G. These peaks 
are separated by twice the Zeeman energy for SIM (c), and are closer for SDM (f). 
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may be essential especially on the normal- metal side, are less important in that regime). 

Another feature differentiating the SIM case from the SDM situation is the absence of 
the spin-up features of conductance spectra for high magnetic fields and for the FF state. 
It is difficult to say, if this feature is model-independent or characteristic to the model with 
dispersion relation of a free-electron gas with renormalized masses. 

Andreev reflection spectroscopy in magnetic field has already been reported in a few 
compounds.—"— For example in MosSby point contact AR spectroscopy lead to identification 
of this compound as an unconventional superconductor.— Such measurements have also 
been performed on pure and Cd-doped CeCoIns.— This compound, as a heavy- fermion 
superconductor and possibly host to the FFLO phase, is a natural candidate for verification 
of the present results. Spectra presented in Fig. 4 of Ref. |75| resemble our Fig. |9^, with 
splitting between the spin-up and the spin-down features of the order of 8 T in fields of 
approximately 2 T. This might indicate that h^or TT h {hcor enhances h), but in CeCoIus 
the one-band model assumed in our calculations may not be sufficient^ and therefore, our 
interpretation is only a speculation.— On the other hand, for a two-band model with strong 
correlations the hcor terms are also present (for both bands), and our conclusions should also 
hold. 

Let us note that, in view of the present results, the AR spectra for the case of the 
BCS state with (100) contact and in the tunneling limit (high Z) would be most helpful 
in detecting the effect of strong correlations in superconductors. Such configuration can 
be studied by both Andreev reflection spectroscopy of a planar junction, as well as by the 
Scanning Tunneling Spectroscopy technique. 



VI. CONCLUSIONS 

In this paper we have provided a detailed analysis of the conductance spectra of a normal 
metal - strongly-correlated superconductor junction. The splitting of conductance peaks in 
the strongly correlated case differs from that in the uncorrelated case. It is equal to twice 
the Zeeman energy only in the latter case and in the correlated case it may be smaller or 
larger depending on the details of the electronic structure. We identify this feature as one of 
the hallmarks of strong correlations in the superconducting phase, as it should hold true for 
other models with different dispersion relations. It is most clearly visible in the case of BCS 
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superconductor with (100) contact and in the tunnehng regime (high Z). In other cases it 
is also present, but the spin-resolved conductances must be analyzed in order to identify the 
splitting unambiguously. 

It would be interesting to examine other spectroscopic methods, such as the Josephson 
tunneling in the SQUID geometry for the systems with strong correlations (and specific 
features resulting from them: the spin-dependent masses and the correlation field). Such 
analysis should be carried out separately as it may lead to a decisively distinct interference 
pattern in an applied magnetic field. 
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